Two-loop self-energy correction in high-Z hydrogen-like ions 
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The calculation of the two-loop self-energy correction to the Lamb shift is at present one 
of the most challenging problems in bound-state QED. Until very recently, this project has 
\ been addressed to mainly within the Za-expansion approach. In it, the two-loop self-energy 

H 



■ contribution is represented as an expansion over Za and ln(Za;), 



F{Za) = + iZa)B5o + {Zaf 
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+L2S62 + L56i+56oJ +■■■ , (1) 

where F{Za) = A£'/[m(a;/7r)^(Za)^], L = \n{Za)~'^. Whereas the lowest-order term B40 
has been known for a long time, calculations of higher-order contributions have not been 
accomplished until recently. The correction B^q was found to be surprisingly large 
B^Q = —24.27. This result has significantly changed the theoretical prediction for the Lamb 
shift in hydrogen and resolved the disagreement with the experimental value existing at 
ihat time. The leading logarithmic term Bq^ was derived first in 0] and later confirmed in 
. (It should be noted that the first evaluation 3], while yielding the right result, is not 
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completely correct, as will be discussed below.) The two remaining logarithmic corrections 
Bq2 and Bqi have also been elaborated lately by Pachucki P]. Again, as in order {Za)^, the 
result obtained turned out to be surprisingly large. The numerical value of Bqi is 50.3, which 
reverses the sign of the overall logarithmic contribution for hydrogen. This indicates that the 
convergence of the Za expansion for the two-loop self-energy correction is remarkably slow, 
and a conclusion has been drawn in |^ that a numerical evaluation with Dirac-Coulomb 
propagators is desirable even for hydrogen. 

The calculation of the two-loop self-energy diagrams (Fig. ^ without an expansion in 
Za started with the irreducible contribution of the diagram (a) (known also as the loop- 
after-loop (LAL) correction), which is by far the simplest part of the total set. Such an 
evaluation was first accomplished in for high-Z ions, and later in [3] for all ions, including 
hydrogen. The latter investigation demonstrated a rather peculiar behaviour of the LAL 
correction in the \ow-Z region. It was shown that for hydrogen its actual value was of a 
different sign and magnitude than the value based on first two terms of the Za expansion. 
In addition, a different result was found in j3| for the leading logarithmic contribution Sga 
as compared to the analytical evaluation . (We note that in the latter work the Bqs term 
was evaluated for the whole set of two-loop self-energy diagrams. However, it was argued 
that it originated from the LAL contribution only.) As a result, a question was raised in 0] 
about the possibility that the Za expansion for the two-loop self-energy could be inadequate 
even for hydrogen. This speculation attracted attention and several investigations followed. 
The subsequent numerical calculation ^ claimed to be compatible with the analytical result. 
However, the third numerical evaluation by one of us joj confirmed the first result . At the 
same time, the total value of the -Bgs contribution was confirmed independently by several 
groups, e.g., in To throw light on this intricate situation we performed an analytic 
calculation of the Bq^ term separately for the LAL correction and found agreement both 
with our previous numerical result and with that of Our conclusion was that the LAL 
correction provided an additional cubed logarithmic contribution that had been omitted in 
the original analytical calculation j^. However, this additional term vanishes when the whole 
set of two-loop self-energy diagrams is taken into account. Recently, analogous additional 
terms were reported for the leading logarithmic contribution for P states 11 1. 

The evaluation of the remaining contributions in Fig. is by far more difficult. These 
contributions are: the reducible part of the diagram (a), the overlapping diagram (b), and 
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the nested diagram (c). The first attempt to evaluate them to all orders in Za was made by 
Mallampalli and Sapirstein In that work, the contribution of interest was rearranged 
in 3 parts, referred to by the authors as the "M", "P", and "F" terms. (We will discuss 
this separation in more detail below.) Mallampalli and Sapirstein calculated only the M 
and F terms, while the P term was left out since a new numerical technique had to be 
developed for its computation. In addition, since the numerical procedure turned out to be 
very time consuming, the actual calculation of the M term was carried out only for two ions, 
uranium and bismuth. Subsequently, in the investigation by two of us j3l we accomplished 
the computation of the remaining P term for Z = 83, 90, and 92, which formally completed 
the calculation of the two-loop self-energy. However, as we will see, the rearrangement of the 
whole correction into the M, P, and F terms is artificial since all the three are divergent. 
A proper treatment should include these terms simultaneously. In addition, more than two 
points (in Z) are needed in order to analyze the Z dependence of the correction and to 
compare it with the known terms of the Za expansion. All these issues are addressed to in 
the present investigation. 

Let us now turn to the evaluation of the two-loop self-energy diagrams. The first problem 
to be solved is the separation of ultraviolet (UV) divergences. The standard method of 
renormalization in QED is developed for diagrams involving only free-electron propagators, 
treating them in momentum space. Therefore, our strategy is to subtract similar diagrams 
with electron propagators containing zero or one interaction with the nuclear Coulomb field 
in order to make the corresponding point-by-point difference UV finite. The subtracted 
diagrams can be then evaluated in momentum space or in the mixed momentum-coordinate 
representation. For the first-order self-energy, this approach was first implemented in jl^ . 
The situation is much more difficult in the case of the two-loop self-energy. Here, for the first 
time, we encounter the overlapping UV divergences. For example, the diagram in Fig. ^b) 
can be considered as consisting of two overlapping vertex subdiagrams, each of which is UV 
divergent. The presence of the overlapping divergences makes the structure of subtraction 
terms much more elaborate than that in the first order. Moreover, some of these terms 
contain both bound-electron propagators and UV-divergent subdiagrams. Such situation 
had never been encountered before, and a new numerical technique had to be developed for 
the evaluation of these subtraction terms. 

Following we rearrange the contribution of the diagrams in Fig.C]in 4 parts: the LAL 
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part, the M, P, and F terms. The LAL correction is defined by the irreducible part of the 
diagram (a). Since its evaluation is relatively easy and has been performed by several groups, 
we do not discuss it here. The M term is diagrammatically represented by Fig. [21 It consists 
of 3 parts originating from the nested diagram, the overlapping diagram, and the reducible 
part of the diagram (a). The subtractions in the M term are chosen so that each of these 3 
parts is separately UV finite. Next, we should account for the subtracted terms. Those that 
contain only free-electron propagators can be treated in momentum space using the standard 
Feynman-parametrization technique. For those that involve bound-electron propagators, we 
introduce additional subtractions that remove the overlapping UV divergences. This is 
graphically represented by Fig. IHl The corresponding contribution is referred to as the P 
term. It consists of 3 parts, each containing only single UV-divergent subgraphs. Taking 
the first part as an example, we see that the difference shown in the figure is UV divergent 
only due to the inner self-energy loop, while the divergence due to the outer self-energy loop 
is canceled. Finally, we collect all terms we have subtracted and denote them as the F term 
depicted in Fig. 0] It consists of Feynman diagrams that contain free-electron propagators 
only. 

We should also mention the infrared (IR) reference-state divergences that are present in 
the M and P terms. These singularities can occur in bound-state QED calculations when 
energies of the intermediate states in the spectral decomposition of electron propagators 
coincide with the valence-state energy. An analysis given, e.g., in shows that the IR- 
divergent terms cancel each other in the sum of the M and the P term. To sum up our 
discussion of divergences, we separately write divergent contributions to the terms under 
consideration: 

AEm = AEIj - AEiR , (2) 
AEp = AE^ + AEiR + lWAE^|+\ (3) 
AEf = AEi + B^'^AE^^\ (4) 

where the index / labels finite contributions, AEi^ is the IR-divergent contribution, L^^^ and 
B^^^ are the one-loop renormalization constants that fulfill the Ward identity L^^'^ = —B^^\ 
and AE^E is the many-potential part of the one-loop self-energy correction. 

We now turn to the numerical evaluation of these terms. It was carried out in the 
Feynman gauge. The P term was evaluated along the lines described in detail in our previous 



investigation The calculation of the F and M terms was performed in a way, in many 
respects similar to that of |12']. The details of the calculation will be published elsewhere. 
Here, we focus on major novel features of our evaluation as compared to Q|. The first point 
is a different treatment of the reference-state IR divergences. In ^2!, they were regulated 
by altering the valence energy Ea to Ea = SaiX — 5)- The actual calculations were performed 
keeping a finite regulator 5, and the limit 5—^0 was evaluated numerically. According to our 
experience, that approach, while being technically easy to handle, does not allow to control 
the accuracy of the computation effectively. In our approach, we introduce some subtractions 
in order to make the terms under consideration finite, separating IR divergences in the form 
Aii^iR. The divergent contributions cancel each other explicitly in the sum, and we can 
perform the whole computation without introducing any actual IR regulator. However, in 
order to allow the term-by-term comparison with the previous evaluation 12], we performed 



our calculations with the regulator 5 as well. 

The second new feature of our approach is a different procedure employed for the double 
summation over the partial waves in the M term. In the photon angular momenta li 
and I2 were chosen as independent expansion parameters. We found it technically more con- 
venient to employ for this purpose the absolute values of the relativistic angular parameter 
K of two electron propagators, |k| = j + 1/2. Thus, we turn the nested and overlapping 
contributions to the M term into tables of values |k2|. Next, we perform a resummation 
of the table: X\t^^\y^\ Yij, where i = \\k.i\ — \k2\\, j = + \i^2\- Finally we sum up the 
table: first, we fix i and extrapolate the sum over j to infinity, and then sum over i and 
estimate the tail of the expansion. 

Now we discuss the computer time necessary for the evaluation of the M term. In the 
previous evaluation by Mallampalli and Sapirstein, a total time of 7323 h was required for 
a given value of Z. In our numerical approach, the typical time of the evaluation of one 
element X^^,^^ is about 1 h for the IBM PWR3 processor with 350 MHz, both for the nested 
and the overlapping diagram. The typical number of elements for a given Z was 440 for the 
nested diagram and 320 for the overlapping diagram. This shows that the time consumption 
in our numerical procedure is smaller than that of although it is still very large. 

In Table I we present the numerical results for finite parts of the M, P, and F terms. 
The table shows that our numerical values for the LAL and F terms agree very well with the 



ones from 



12| | but there is a significant deviation for the M term. More specifically, our 
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calculation for Z = 92 yields —2.137, 4.679(2), and —3.837(2) for the reducible, nested, and 

overlapping contributions to the M term, respectively. These results should be compared 

I I 

correspondingly with —2.137, 4.669(5), and —4.387(5) from 12]. We see that the leading 
source of discrepancy is the overlapping diagram. Taking into account the complexity of the 
computation, it is difficult to suppose what the reason for this disagreement could be. 

As in the case of the one-loop self-energy, the evaluation becomes problematic very fast as 
Z decreases. It is due in part to the fact that some individual contributions exhibit a nearly 
Z-independent behaviour, while the total correction scales as (Za)^. Numerical problems 
restricted our calculation to the region Z > 40. In Fig. El we compare our non-perturbative 
results with the known terms of the Za expansion. Although we can not as yet say anything 
conclusive about the higher-order terms, the figure suggests that the results obtained by two 
different methods could be compatible. 

Summing up, we have evaluated all contributions to the two-loop self-energy correction for 
H-like ions with Z > 40. As this correction has been the major source of the uncertainty of 
theoretical values for the ground-state Lamb shift in these systems, our evaluation improves 
their accuracy by an order of magnitude Q]. While the experimental precision for H- 
like uranium is not presently sufficient to probe the new contribution, this should become 
possible when the experiment currently planned at GSI jlS^j is completed. The question of 
utmost importance is to extend the present evaluation to low-Z ions, where higher-order 
terms could enter at the level of experimental interest even at Z =1 j^, as well as to excited 
states. For the 2pi/2-2s transition in Li-like high-Z ions the two-loop self-energy presently 
defines the uncertainty of the theoretical prediction and can be probed by comparing 
with experimental data available. 
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FIG. 1: Two- loop self-energy diagrams. Double line indicates an electron propagating in the 
Coulomb nuclear field. It is understood that the corresponding mass counterterms are subtracted 
from the diagrams. 
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FIG. 2: Diagrammatic representation of the M term. The dashed line denotes the interaction with 
the Coulomb nuclear field. AE'gE is the first-order self-energy correction. 
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FIG. 3: Diagrammatic representation of the P term. The last part should be counted twice, 
accounting for two equivalent terms. 
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FIG. 4: Diagrammatic representation of the F term. The last diagram on the right in the first two 
rows should be counted twice, accounting for two equivalent diagrams. 
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FIG. 5: The results of our numerical evaluation to all orders to Za (dots) together with the first 
two terms of the Za expansion (solid line) and all known terms of the Za expansion (dashed line). 
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TABLE I: Individual contributions to the two-loop self-energy correction expressed in terms of 
F{Za). 



z 


LAL 


F term 


P term 


M term 


Total 


40 


-0.871 


19.50 


-30.13(15) 


10.50(18) 


-1.00(26) 


50 


-0.973 


10.03 


-14.42(7) 


4.04(7) 


-1.33(10) 






10.02" 








60 


-1.082 


5.72 


-7.48(4) 


1.21(2) 


-1.63(4) 


70 


-1.216 


3.497 


-4.03(3) 


-0.14(1) 


-1.89(3) 




-1.216^ 










83 


-1.466 


1.938 


-1.831(13) 


-0.990(5) 


-2.349(14) 






1.937'^ 




-1.66(1)" 




92 


-1.734 


1.276 


-1.030(9) 


-1.295(3) 


-2.784(10) 




-1.733^ 


1.274" 




-1.855(7)" 




100 


-2.099 


0.825 


-0.635(6) 


-1.473(3) 


-3.382(7) 






0.825 


n 


n 
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